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ABSTRACT 

We investigate the radial dependence of the spectroscopic properties, in particular the Mg2, 
<Fe> and H/3 spectroscopic indices, in a sample of galaxies spanning a wide range of absolute 
luminosity in the Coma cluster. After allowing for the magnitude dependence of these indices, 
we find a significant gradient in Mg2, in the sense that galaxies in the core of the cluster have 
stronger Mg2. We find only weak gradients in <Fe> and H/3. Using the model grids presented in 
an earlier paper in this series, we attribute the Mg2 gradient to changes in metal abundance. One 
possible mechanism to create this abundance gradient is pressure confinement by the intracluster 
medium of material from Supernova driven winds early in the history of the galaxies. 



Subject headings: galaxies: clusters — galaxies: clusters: individual(Coma) 
lenticular — galaxies: evolution 



galaxies:elliptical and 



1. Introduction 

This is the fifth paper in a series which exam- 
ines the properties of galaxies in the Coma cluster 
in the magnitude range 12 < R < 20 and in a wide 
variety of environments. In this paper we exam- 
ine the spatial dependence of line indices which 
depend upon metal abundance and age, and of 
the luminosity weighted mean stellar abundances 
and ages which we derived from these indices in 
Paper III of this series (Poggianti et al. 2001). In 
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this paper we examine the spatial dependence of 
the line indices diagnostic of metal abundance and 
age and investigate the origin and implication of 
such gradients. 

The properties of galaxies in clusters have long 
been known to depend strongly upon the environ- 
ment. This was shown most clearly by the survey 
of Dressier (1980a, b) who showed a dependence of 
galaxy morphology upon environment in the sense 
that the proportion of E and SO galaxies increases 
with increasing density. This has been taken as 
evidence for a merger process occurring in clus- 
ters which transforms disk galaxies to spheroids 
at an early epoch. Further work on clusters at 
higher redshifts (Dressier et al. 1997) shows that 
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this dependence is qualitatively similar in highly 
concentrated clusters at earlier epochs. 

The importance of the intracluster medium has 
become clear more recently, the hot gas mass in 
clusters is typically 5 times the luminous mass 
in galaxies, and that gas typically has a metal 
abundance (specifically iron abundance) of 0.2 - 
0.5 solar (Finoguenov et al. 2000; De Grandi & 
Molendi 2000). The majority of the heavy ele- 
ments in a cluster are in the hot gas, but all of 
these heavy elements must have come from galax- 
ies at some stage. Finoguenov & Ponman (1999) 
and Finoguenov et al (2000) study the distribu- 
tion of both iron and a elements in the hot gas 
component in clusters, and find that there are ra- 
dial gradients in the abundances, stronger in iron 
than in a elements. They conclude that Type 
II supernovae dominate the enrichment process in 
the outer regions of clusters, but in their centers 
Type I supernovae are more important. However 
in the Coma cluster, which is the subject of this 
study, Arnaud et al. (2001) find that the abun- 
dance profile is flat at 0.25 solar out to 15 arcmin- 
utes radius, except for a region of enhanced abun- 
dance less than 30 arcsec in radius ccnterd on NGC 
4874. These authors present only a single abun- 
dance profile, presumably derived assuming solar 
abundance ratios. 

As the galaxies are the source of the heavy 
elements, it is important to study and under- 
stand their properties. Line indices and luminos- 
ity weighted ages and abundances for samples of 
Coma cluster galaxies have been published by Jor- 
gensen (1999) and Terlevich et al. (2000), but our 
sample has a much greater range in luminosity, 
and covers a wide range of radius and of environ- 
ment. 

2. Observations 

Our data are taken from photometric and spec- 
troscopic observations of two 40 x 30 arcminute 
fields in the Coma cluster, with the Tokyo Mosaic 
CCD camera and the WYFFOS fibre fed spectro- 
graph respectively. Details of the observing setup 
and data reduction arc given in earlier papers in 
this series. Photometric observations and reduc- 
tion are described in paper I (Komiyama et al. 
2001). Paper II (Mobashcr ct al. 2001) describes 
the spectroscopic sample selection. In paper III 



(Poggianti ct al. 2001) we derive line indices on 
the Lick system and luminosity weighted ages and 
metallicities from these indices. In the current pa- 
per we analyse the following quantities: galaxy 
position, R magnitude, (B-R) color and concen- 
tration index from Paper I; and the spectroscopic 
indices Mg2, <Fe> and 11/3 from paper III, to- 
gether with luminosity weighted ages and metal- 
licities derived from the the model grids presented 
in the Mg2 - H/3 diagram in Figure 4 of that paper. 

Our galaxy sample selection is described in pa- 
per II. The spectroscopic sample comprises two 

subsamples, a bright sample of galaxies at R<18 
with known redshifts from CoUess & Dunn (1996) 
and Colless (1998) and a faint sample selected 
from our own photometric imaging survey (Papc;r 
I) to have 18<R<20 and 1<B-R<2. As described 
in Paper II, the samples are both believed to be 
free of morphological bias. Our faint spectroscopic 
sample is not complete, but the galaxies observed 
are drawn at random from the photometrically de- 
fined sample. 

The sample discussed here comprises 229 galax- 
ies from our spectroscopic survey for which the 
formal errors on the Mg2, <Fe> and H/3 indices 
are less than 0.04 magnitude, 1.5 A and 1.5 A re- 
spectively, and which did not show emission lines 
on visual inspection. Of these, 163 are in the cen- 
tral Coma 1 field and 66 in the Coma 3 field which 
is located to the southwest of the cluster center, 
described in Paper II. If more than one such spec- 
trum was available then the line index values were 
combined with weighting proportional to the in- 
verse of the variance on each line index. Because 
we require small errors in the indic(!s for inclu- 
sion in the sample we discriminate against faint 
galaxies and in particular low surface brightness 
galaxies, for low surface brightness the signal-to- 
noise ratio of our spectra is often too poor for in- 
clusion. Rejecting emission lines (which otherwise 
will dilute the Balmer indices in particular) dis- 
criminates against galaxies with ongoing star for- 
mation. 

In this paper we examine the dependence of 
galaxy properties on the distance of the galaxy 
from the the cluster center, which has been taken 
to be at RA = 12h 59m 42.8s, Dec = 27° 58' 14", 
Equinox J2000.0, a position near the peak of the 
X-ray emission (White et al. 1993) and some 1.8 
arcminutes East of NGC 4874. However we do 
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not completely sample the azimuthal distribution 
at radii greater than 15 arcminutes. Our outer 
field is to the South West of the center, and in- 
cludes NGC 4839 and the group of galaxies aroimd 
it, which are believed to be falling into the Coma 
cluster (e.g. Neumann et al. 2001). The effect of 
the inclusion of this group in our sample is consid- 
ered later in the paper. 

3. Radial dependence of line indices 

In Figures 1 - 3 we plot the Lick Indices Mg2, 
<Fe> and H/9 as a function of radius in the cluster. 

The Mg2 index (Figure 1) shows a clear trend 
with radius, galaxies in the outer part of the clus- 
ter have lower Mg2 indices than those in the core. 
There is a substantial scatter at all radii which 
is due to the strong trend of this index with 
galaxy magnitude (paper III); the fainter galax- 
ies have weaker Mg2, and larger errors, because 
of the lower signal-to-noise of the spectra. The 
lines show weighted (dashed line) and unweighted 
(dot-dash line) least squares fits to the data, there 
is a large offset between the lines because the 
brighter galaxies which have stronger Mg2 have 
better signal-to-noise and hence lower errors and 
higher weights in the weighted fit. However both 
fits show the trend which is visible in the data; to- 
wards lower Mg2 at larger radii within the cluster. 

In Table 1 we present the statistical properties 
of the fits, to show the significance of the correla- 
tions. The important statistics in this context are 
"Student's" t, and P, the probability of obtaining 
a value o{\t\ this high by chance in an uncor related 
sample. 

The equivalent plots for the <Fe> and H/3 in- 
dices are shown in Figures 2 and 3 respectively, 
the <Fe> index shows a weaker gradient than 
Mg2, perhaps masked by higher errors, whereas 
H/J shows a marginal positive gradient. 

Although Figure 3 shows marginal evidence for 
a positive gradient in 11/3 the statistical parame- 
ters in Table 1 show that this is not significant. In 
this Figure we do not see evidence for a concen- 
tration of Balmer line strong galaxies in the NGC 
4839 region of the Coma cluster, as was proposed 
by Caldwell ct al. (1993), although our survey 
area does cover this region. NGC 4839 is 41.7 ar- 
cminutes from our defined cluster center, and we 
see no concentration of H/3 strong galaxies at this 
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Fig. 1. - Mg2 index against radius for our sample 
of galaxies in the Coma cluster. The dashed and 
dot-dash lines are weighted and unweighted least 
squares straight line fits to the data respectively. 



<Fe> agoinst Radius 




Fig. 2. — <Fe> index against radius for our sam- 
ple of galaxies in the Coma cluster. The dashed 
and dot-dash lines are weighted and unweighted 
least squares straight line fits to the data respec- 
tively. 



3 



Fig 


Ordinate 


Abscissa 
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Intercept 
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Mg2 




Radius 


229 


-0.22 


-6.51x10"* 


1.91x10" 


-4 


0.229 


6.03x10-3 


-3.41 


0.08 
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<Fe> 




Radius 


229 


-0.08 


-2.14x10"^ 


1.69x10" 


-3 


2.627 


0.054 


-1.26 


20.7 


3 


H/3 




Radius 


229 


0.19 


6.05x10-3 


2.03x10" 


-3 


1.649 


0.064 


2.99 


0.31 


4a 


Mg2 




R mag 


229 


-0.74 


-3.58x10-2 


2.15x10" 


-3 


0.7385 


3.16x10-2 


-16.7 


< 10"3 


4b 






Radius 


229 


-0.25 


-5.01x10"* 


1.27x10" 


-4 


1.21x10-2 


4.00x10-3 


-3.95 


0.01 


5 


(5(B-R) 




Radius 


229 


-0.45 


-4.22x10"'^ 


5.63x10" 


-4 


0.1052 


1.82x10-2 


-7.51 


< 10"3 


6a 


<Fe> 




R mag 


229 


-0.56 


-0.243 


0.024 




6.149 


0.348 


-10.3 


< 10"3 


6b 


<5(<Fc>) 




Radius 


229 


-0.07 


-1.37x10-3 


1.40x10" 


-3 


3.34x10-2 


4.45x10-2 


-0.98 


33 


7a 


up 




R mag 


229 


0.31 


0.159 


0.032 




-0.541 


0.475 


4.94 


< 10"3 


7b 


s(m 




Radius 


229 


0.19 


5.45xl0"3 


1.93x10" 


-3 


-0.133 


0.061 


2.82 


0.5 


8 


Metallicity 




Radius 


208 


-0.16 


-4.34xl0"3 


1.89x10" 


-3 


-0.5028 


0.0603 


-2.30 


2.2 
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Metallicity 




R mag 


208 


-0.58 


-0.1967 


0.0193 




2.619 


0.319 


-10.2 


< 10"3 


9b 


(5(Metallicity) 




Radius 


208 


-0.11 


-2.46xl0"3 


1.55x10" 


-3 


6.06x10-2 


4.95xl0"2 


-1.59 


11.4 


11 


c5(Metallicity) 




Radius 


143 


-0.22 


-4.20xl0"3 


1.58x10" 


-3 


9.39xl0"2 


4.55xl0"2 


-2.66 


8.8 



Table 1: Table showing the parameters describing the fits illustrated in Figures 1-9 and 11. Column 1 
refers to the figure in the paper, columns 2 and 3 list the dependent and independent variable respectively. 
Column 4 gives the number of points used in each fit, and column 5 the correlation coefficient between the 
variables. Column 6 gives the slope of the fitted straight line (illustrated in the plot) and column 7 the formal 
error on this value. Columns 8 and 9 are the intercept with the vertical axis, and its error, respectively. 

Column 10 gives "Student's" t statistic, which is related to the correlation coefficient by t = i^^i i and 
is also identically equal to the ratio of the slope and its error. Column 11 gives the two-tailed (significance) 
probability, expressed as a percentage, of a value of \t\ this high occurring by chance. 
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Fig. 3. — H 13 index against radius for our sample 

of galaxies in the Coma cluster. The dashed and 
dot-dash lines are weighted and unweighted least 
squares straight line fits to the data respectively. 



radius. Caldwell ct al. measured US rather than 
H/3, which is less affected by emission. Emission 
line galaxies have been eliminated from our sam- 
ple following visual inspection, although some low 
level of emission may remain (possibly giving rise 
to the negative outliers in Figure 3). Again, these 
galaxies appear to be distributed evenly through- 
out the cluster. 

3.1. Correcting for the luminosity depen- 
dence of the indices 

In Paper 3 we show a tight correlation between 
the metallicity derived from the Mg2 index and 
galaxy luminosity for our sample. Here we investi- 
gate whether this dependence alone is sufficient to 
explain the radial gradient that we see in Figure 1. 
In Figure 4a we plot Mg2 against R magnitude for 
the sample used here. This shows the well defined 
correlation between Mg2 and magnitude, which is 
discussed in Paper 3. To take this relationship out 
of the radial dependence, in Figure 4b we plot the 
residuals in the Mg2 index from a weighted least 
squares straight line fit against radius. 

The correlation with radius is clearly still 
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Fig. 4. — a) Mg2 against R magnitude for our sample, the dashed hne is a weighted least squares straight 
hne fit to the data, b) Mg2 residual from the fit in (a) against radius within the cluster. The dashed line is 
a weighted least squares straight line fit. 



present in the residuals, and the scatter has been 
reduced, which shows that the correlation in Fig- 
ure 1 is not entirely due to index-magnitude rela- 
tion. 

A similar effect is seen if we plot the residu- 
als from the Color-Magnitude relation (paper I) 
against radius (Figure 5). 

In Figures 6 - 7 we repeat this procedure for 
our mcasurcmcints of the <Fc> and H/3 indices, to 
see whether we can identify residual radial trends 
in these indices which may have been masked by 
the scatter introduced by the index-magnitude re- 
lations. <Fe> shows a strong dependence upon 
magnitude (Paper III and Figure 6a), and when 
this relation is taken out we sec no dependence of 
the residual upon radius (Figure 6b). 

This result is surprising given the strong gradi- 
ent seen in Mg2, and if confirmed suggests a dif- 
ferent radial dependence of the enrichment mecha- 
nisms for these two elements. However we have to 
consider the possibility that it is simply the larger 
scatter in <Fe> which masks a radial trend. For- 
mally the slope in Figure 4b is -0.00050 ± 0.00013 
mag arcmin"^ and that in Figure 6b is -0.00137 
± 0.00140 A arcmin . If the slopes of the resid- 
ual/radius relations were to be in the same ratio 
of the slopes of the index/magnitude relations for 
Mg2 and <Fe> (Figures 4a and 6a) then we would 



Colour-Residual against radius 




20 40 60 80 

Rodius (arcmin) 



Fig. 5. — Residual from the color-magnitude re- 
lation for the galaxies in the spectroscopic sam- 
ple, plotted against radius within the cluster. The 
dashed line is a least squares straight line fit to 
the data. 
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Fig. 6. — a) <Fe> against R magnitude for our sample. The dashed line is a weighted least squares straight 
line fit to the data, b) <Fe> residual from the fit in (a) against radius within the cluster. The dashed line 
is a weighted least squares straight line fit. 



expect a slope in Figure 6b of -0.0034 A arcmin ^ , 
which is 1.4(7 away from the observed value. 

There is a very small positive radial gradient in 
the H/3 residual, amounting to only 0.4 A over the 
entire radial range. 

3.2. Metallicity or age as a driver of the 
radial dependence of the indices 

In paper 3 we interpreted our measurements of 

Lick indices in terms of luminosity weighted stellar 
age and metallicity with the help of Single Stellar 
Population models described in that paper. The 
models show that the Mg2 index in particular is 
sensitive not only to the metallicity of the popu- 
lation, but that at high metallicity it is also sensi- 
tive to the age of the population. This problem is 
sometimes referred to as the "age-metallicity de- 
generacy" . However, if galaxies really did consist 
of a single age and metallicity population, then a 
pair of indices such as Mg2 and H/3 would yield an 
unambiguous estimate of the age and metallicity 
of the population, subject to observational error 
and the correctness of a particular set of models 
such as those we use in Paper 3. 

The problems with interpretation of the Lick 
indices are largely due to the uncertainties in the 
models. First, the models are only Single Stel- 
lar Population models, and a galaxy such as our 



own does not consist of stars with a single age and 
abundance. Second, our understanding of stellar 
evolution is as yet incomplete. The limitations 
of current models are most apparent from our 
lack of understanding of Galactic globular clus- 
ters. We still don't know if age is the dominant 
second parameter affecting horizontal branch mor- 
phology (e.g. Chaboyer et al. 1996; Stetson et 
al 1999; Rich et al. 1997), and some authors 
claim that the age differences between clusters of 
similar metallicity are too small to explain their 
difference in HB morphology. A related prob- 
lem is the embarrassing discrepancy between the 
spectroscopic and CMD ages for 47 Tuc (Gibson 
et al. 1999), though population synthesis mod- 
els with improved physics reduce the discrepancy 
(Vazdekis et al. 2001) or possibly eliminate it 
completely (Schiavon et al. 2001). Certainly, if 
we cannot model the relatively simple globular 
clusters, which are coeval and chemically homo- 
geneous, we will have little hope of understanding 
more complicated stellar systems such as galaxies. 

Accepting these caveats, we use the models de- 
scribed in Paper III to investigate whether the ra- 
dial gradients in the indices we observed are due 
to metallicity or age gradients, or a combination 
of the two. 

In Figure 8 we plot the metallicity from the 
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Fig. 7. — a) H/? against R magnitude for our sample, 
line fit to the data, b) H/3 residual from the fit in (a) 
a weighted least squares straight line fit. 

grids in the H/3/Mg2 diagram in Paper 3 against 
radius, there is a gradient of around 0.4 dex from 
the center of the cluster to 80 arcminutcs radius. 

In Figures 9a and 9b we again take out the de- 
pendence of metallicity upon luminosity, Figure 
9a is the correlation between metallicity and R 
magnitude, and in Figure 9b we plot the residuals 
from the fit to this relation against radius. When 
the strong luminosity dependence is taken out we 
find a gradient Alog[Fe/H] ~ -0.2 between the cen- 
ter of the cluster and 80 arcminutes radius. The 
significance of this gradient is not as well estab- 
lished as those in the raw indices, Table 1 shows 
an 11% probability of a gradient of this size aris- 
ing by chance, i.e. it is established at the 89% 
confidence level. Further observations of a larger 
sample of galaxies are required to firm up this re- 
sult. 

A feature of Figure 9a is the lack of galaxies 
in the lower left quadrant of this plot, i.e. there 
arc very few bright metal-poor galaxies, although 
there are rather more faint metal-rich ones, al- 
though our spectroscopic sample at R < 16 is com- 
plete. This distribution, together with the very 
strong correlation observed in this diagram, sug- 
gests that processes which form smaller galaxies 
from larger ones (tidal stripping, fragmentation 
and harassment) were more important in the evo- 
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The dashed line is a weighted least squares straight 
, against radius within the cluster. The dashed line is 



lution of the cluster than processes which form 
bright galax;ies out of faint ones (e.g. near equal- 
mass mergers), at least at times since the bulk of 
the stars were formed. 

In Figure 10 we plot histograms of the derived 
luminosity weighted stellar age from the models 
in Paper III, divided into samples at projected 
radius greater than and less than 40 arcminutes. 
Our ages are derived from the right hand panel 
of Figure 5 of that paper, as pointed out there, 
many galaxies lie below the grid lines, i.e. at lower 
H/3 strength than any models predict, these galax- 
ies are allocated an age of 30 Gyr. The cause of 
this discrepancy is discussed in Paper III. Figure 
10 shows that a higher proportion of the galax- 
ies in the inner sample than in the outer lie be- 
low the grid. The ages attributed to these par- 
ticular galaxies are not satisfactory for address- 
ing the problem addressed in this paper, and be- 
cause the significance of the difference in the age 
distributions outside and inside 40 arcminutes de- 
pends critically upon our treatment of these galax- 
ies, we can only test such differences using non- 
parametric rank tests. Using the two-sample form 
of the Kolmogorov-Smirnov test, we test the sam- 
ples against the null hypothesis that they are 
drawn from the same distribution. Taking at face 
value the rank of the ages of the galaxies in Fig- 
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Fig. 9. — a) IVIetallicity estimated from the model grid in the H/3/Mg2 diagram from Paper 3, plotted against 
R magnitude. The dashed line is an unweighted least squares fit. b) Residuals from the fit in (a) plotted 
against radius within the cluster. 



ure 10, and using the two-sample form of the test, 
we find that the distributions differ at the 93% 
confidence level. However ignoring those galax- 
ies placed at an age of 30 Gyr, we find that the 
age distributions for those galaxies with well es- 
tablished ages differ only at the 30% confidence 
level, i.e. those distributions are entirely consis- 
tent with each other. 

Maraston & Thomas (2000) present composite 
population models using isochrones which repro- 
duce the strong Balmer lines seen in old, metal- 
poor galactic globular clusters. Their models of 
early type galaxies consist of two old populations, 
one metal poor and one metal rich. They propose 
that these models can give rise to the observed 
scatter in cluster galaxy H/3 strengths as a combi- 
nation of different proportions of old metal poor 
and old metal rich populations without the need 
for a young population. If there were a metal poor 
population in the galaxies in our sample which be- 
came more important further out in the cluster, 
then this could mask an age gradient, but this age 
gradient would be in the sense that galaxies at 
the center of the cluster were younger. There is 
no firm evidence for such an effect, and although 
the work of Maraston & Thomas (2000) suggests 
that there may be galaxies in our sample for which 
we are deriving unrealistically young ages, these 



galaxies seem to be spread evenly throughout the 
cluster. 

3.3. Dependence of the abundance gradi- 
ent on the morphology and age 

It is well known that there is a dependence of 
galaxy morphology upon environment within clus- 
ters (Dressier 1980a, b), so in this section we con- 
sider whether the mctallicity gradient which we 
see is the same effect. At first sight this seems un- 
likely, as the processes which give rise to transfor- 
mations between morphological types (ram pres- 
sure stripping, galaxy harassment) are unlikely to 
affect the galaxy metallicity, and any effects on 
the mctallicity-magnitude relation wiU result from 
effects on the luminosity, and is second order to 
effects upon the luminosity function. Moreover 
the processes which determine the galaxy metal- 
licity are either internal and depend primarily 
upon the galaxy mass (resulting in the metallicity- 
magnitude relation), or they are environmental 
and involve transfer of gas between the galaxies 
and the intracluster medium, and are unlikely to 
depend upon morphology. 

Our sample consists largely of galaxies fainter 
than V=17, for which we have no objective or 
subjective morphological classification. Inevitably 
there will be some spiral and irregular galaxies in 
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Fig. 8. — Metallicity estimated from the model 
grid in the H/?/Mg2 diagram from Paper 3, plotted 
against radius within the cluster. As there arc now 
no formal errors on these values the dashed line 
represents an unweighted least squares fit. 



the sample. Tcrlcvich ct al. (2001) plot the (U- 
V) color-magnitude relation for galaxies of vari- 
ous morphological types in the Coma cluster (their 
Figure 6). The color - magnitude (CM) relation 
is driven mostly by the metallicity-mass relation 
(Faber 1973); Terlevich et al. find that most 
spiral, irregular and unclassified galaxies in their 
sample lie on the same CM relation as the E and SO 
galaxies, although there are a few outliers on the 
blue side of the ridge line. However these outliers 
are very blue, around 1 magnitude bluer in (U-V) 
than the ridge line, which is likely to be due to 
ongoing star formation. Such galaxies would have 
emission lines, and would have been rejected from 
our sample. 

Figure 11 shows the residuals in metallicity 
from the metallicity magnitude relation, plotted 
only for those galaxies with a well established 
mean stellar age greater than 6Gyr. There are 
only 143 galaxies in this plot, due to the uncer- 
tainties in the ages. However the correlation with 
radius is still clear, indeed it is stronger and more 
significant than in the complete sample (Figure 
9b); Table 1, so the gradient is stronger in older 



Fig. 10. — Histogram of the luminosity weighted 
stellar age, for galaxies at radii greater than 40 ar- 
cminutes (upper panel) and less than 40 arcmin- 
utes (lower panel). Galaxies which fall below the 
grid (i.e. at lower H/3 strength than any models 
predict) in Figure 5 of Paper HI are allocated an 
age of 30 Gyr, which gives rise to the peaks at this 
age. 

galaxies, or at least galaxies with older stars. 

4. Discussion 

There is clear evidence for a metallicity gradient 
for the galaxies in the Coma cluster, this presents 
itself as a difference of -0.03 magnitude in the Mg2 
index (Fig 4b) or -0.2 in log[Fc/H] (Fig 9b) be- 
tween the value at the center of the cluster and 
that at 80 arcminutes radius. 

There are good reasons to expect that there 
might be a gradient in the metallicity of galax;- 
ies within clusters. Although galaxies will have 
undergone several orbits or crossings within the 
cluster they will not be well mixed, so galaxies 
which form in the outer regions of a cluster will 
spend most of their time there. The metallicity- 
magnitude relation can be understood in terms of 
the loss of processed material in supernova driven 
outflows (Larson 1974, Vader 1984). This loss of 
the heavy elements is counteracted by the grav- 
ity of the galaxy, hence the metallicity-mass and 
metallicity-magnitude relation. But it can also 
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Fig. 11. — Residuals from fit to the metallicity- 
magnitude diagram, plotted against radius within 
the cluster, only for those galaxies with mean stel- 
lar age greater than 6Gyr. 

be counteracted by external pressure of the intra- 
clustcr medium (Silk ct al 1987). Galaxies formed 
in the densest parts of the cluster would be ex- 
pected to retain more of their heavy elements. 
We only detect a significant gradient in Mg2, not 
<Fe>, so it is possible that the mechanism which 
constrained the processed material was more ef- 
fective at a time when Supcrnovae Type II were 
the dominant source of processed material, early 
in the history of the cluster. This is also supported 
by the evidence of Figure 1 1 , that the gradient is 
stronger in galaxies with an older stellar compo- 
nent. 

This picture is complicated by the possible 
modification of the cluster membership by accre- 
tion of subclusters. There is some evidence from 
the X-ray morphology (Neumann et al. 2001) that 
the NGC 4839 group is being accreted by the clus- 
ter. Any such accretion would be expected to di- 
lute any gradient existing in the cluster. 

Guzman ct al. (1992), investigating the (T-Mg2 
relation amongst the bright ellipticals in Coma, 
find an offset between a core sample and a halo 
sample in the sense that the halo ellipticals have 
weaker Mg2 for a given a. These authors, having 



only one index, cannot distinguish between metal- 
licity and age as an origin of this offset, they show 
that age is a possible origin. Rose et al. (1994), 
using different indices, conclude that there is a dif- 
ference in stellar populations between galaxies in 
rich clusters and those in a field environment, in 
that the luminosity weighted mean stellar popula- 
tion is younger and more metal rich in low density 
environments. They attribute this difference to 
an additional intermediate age population in the 
low density environments. They do not investi- 
gate the difference between high- and low-density 
regions in the same cluster. We note that their 
dataset is heterogeneous, particularly in the de- 
tectors used to gather it, and includes not only 
CCD spectra but also some from older image tube 
based detectors. It is not clear whether the dif- 
ference between their conclusions and ours results 
from different observational techniques, or differ- 
ent sample selection, or because the comparison 
between high- and low-density clusters is different 
from the comparison between regions of different 
density in the same cluster. 

Terlevich et al. (2001) find a gradient in the (U- 
V) color residual from the CM diagram with ra- 
dius, in the sense that galaxies in the outer part of 
the Coma cluster are bluer than those in the cen- 
ter. At 30 arcminutcs radius, the galaxies are 0.03 
magnitudes bluer in (U-V) than in the cluster cen- 
ter. Terlevich et al. (2001) attribute this gradient 
to a gradient in the mean stellar age of the galax- 
ies, however (U-V) is also sensitive to metallic- 
ity (Wallerstein 1962; Faber 1973). The relation- 
ship between (U-V) and metallicity in single stellar 
population models is complicated (e.g. Vazdekis 
et al. 1996), and depends upon other parame- 
ters, and indeed differs by worrying amounts be- 
tween models. Inspection of Figure 3 of Vazdekis 
et al. (1996) shows however that a gradient of -0.1 
in log[Fe/II] between and 30 arcminutes radius 
would be sufficient to produce the gradient in (U- 
V) seen by Terlevich et al. without the need for 
an age gradient. 

Balogh et. al (2000) propose a model in which 
star forming galaxies fall into rich clusters at a 
fairly uniform rate over time, and their star for- 
mation is gradually quenched by ram pressure and 
tidal stripping removing their gas over a period 
of a few Gyr. Their model explains the observed 
gradients in star formation rate, as measured by 
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spectroscopic indices in the CNOC cluster sample 
(Balogh et al. 1999). This mechanism would be 
expected to lead to an age gradient in our sample, 
indeed given the small sample and the uncertain- 
ties in the ages of a substantial number of galax- 
ies in Figure 10 we cannot rule this out. How- 
ever as galaxies further out in the cluster would 
have had their star formation truncated on aver- 
age later than galaxies in the core, any effect on 
the metallicity gradient would be in the opposite 
sense to that which we see. We conclude that the 
metallicity gradient is not due to infall of galaxies. 

5. Conclusions 

We detect significant evidence for a negative ra- 
dial gradient in the metallicity of the galaxies in 
the Coma cluster, as measured from pairs of in- 
dices including Mg2. Our <Fe> index measure- 
ments are less accurate, however the gradient that 
we find in this index is significantly shallower that 
we would deduce from the size of the Mg2 gradient. 
Evidence for a gradient in luminosity weighted 
stcillar age is unclear, and we suggest that the gra- 
dient in (U-V) observed by Terlevich ct al. (2001) 
and attributed by them to an age gradient could 
instead be driven by the metallicity gradient that 
we observe. 

Our observations are most easily explained by 
models in which the pressure of the intracluster 
medium confines enriched material in supernova 
driven winds early in the history of the cluster. We 
find no firm evidence for a gradient in the mean 
stellar ages of galaxies with radius. Thus our ob- 
servations support models in which the epoch of 
formation of the bulk of the stars in galaxies does 
not depend upon their location within the cluster. 
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